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When a two-component Bose-Einstein condensate is placed into rotation, a lattice of vortices and 
cores appear. The geometry of this lattice (triangular or square) varies according to the rotational 
value and the intercomponent coupling strengths. In this paper, assuming a Thomas-Fermi regime, 
we derive a point energy which allows us to determine for which values of the parameters, the lattice 
goes from triangular to square. It turns out that the separating curve in the phase diagram agrees 
fully with the complete numerical simulations of the Gross-Pitaevskii equations. We also derive a 
formula for the critical velocity of appearance of the first vortex and prove that the first vortex 
always appears first in the component with largest support in the case of two disks, and give a 
criterion in the case of disk and annulus. 

PACS numbers: 



INTRODUCTION 

When a two-component condensate is set into rotation, 
topological defects of both order parameters are created, 
which lead to more exotic defects than in a single com- 
ponent condensate. Experiments on two component con- 
densates have shown how the condensates can exhibit 
either triangular or square vortex lattices [1]. Accord- 
ing to the values of the interaction strengths, the defect 
patterns can vary a lot, as illustrated in the numerical 
simulations [2, 3]. One specific feature is the appear- 
ance of coreless vortices: the existence of a vortex in 
componcnt-1 corresponds to a peak in component-2 and 
vice-versa. The interaction between vortices and peaks 
leads to changes in the geometry of the vortex lattice. We 
are interested in determining the equations governing this 
vortex peak behaviour in the Thomas-Fermi regime and 
estimating the interaction energy between the lattices of 
the two components. Indeed, for a single condensate, 
the vortex lattice is triangular, while for a two compo- 
nent, the vortex-peak interaction can lead to a square 
lattice. In [4], an asymptotic interaction between two 
half-quantized vortices is derived for two-component ho- 
mogeneous condensates. In this paper, we want to take 
into account the nonhomogeneity of the condensate due 
to the trapping potential and estimate the vortex-peak 
energy according to the parameters of the system. We 
derive an energy depending on the location of vortices 
and peaks and determine for which values of the exper- 
imental parameters, the lattice goes from triangular to 
square. These critical values agree well with the ones 
found from the numerical computations of the full Gross- 
Pitaevskii equations of [3]. We note that in the rapid 
rotation regime, using the lowest Landau level approxi- 
mation, several papers [5, 6] (see [2, 7] for a review) have 
analyzed the transition between triangular to square lat- 



tices. We point out that [8] found a point energy with an 

I I ^ 

interaction term e^'^'^''^' , which is different from ours 
and from [4]. We first review relevant results for a sin- 
gle condensate (see [9]), before moving to the derivation 
of homogeneous equations and the computation of the 
interaction term in two-component condensates. 

For a single component condensate, the wave function 
minimizes the energy 
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= j ^\y^-inxri,\-' + ^{V{r)~n'r'M' + ^\i:\' 

(1) 

under / |V'P = 1, where i7 = flcz is the rotation, V{r) 
is the trapping potential and in most cases V{r) = r^. 
We will denote by Vfj the operator V — zfi x r. For g 
large, at 17 = 0, the ground state 77 of Eg^ approaches 
the inverted parabola 

-(A-r^) 

in the disk of radius = X = gj-K, and goes to 
outside the disk. If the problem is rescaled on a disk of 
size 1, then the analysis of the vortex cores leads to a 
vortex of size and, close to the core, the wave func- 

tion behaves like /(r)e** where / is the solution tending 
to 1 at infinity of 
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/' / 



+ /(l-/2)=0. 



(2) 



This is the equation of a vortex core in a uniform system. 
In the case of a single condensate, from the equation of 
the vortex core, one can estimate the energy of vortices, 
the critical velocity for the nucleation of the first vortex 
and the interaction energy between vortices [9-12] which 
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is 

^Y.\og\p,-p,\+j2\p,\' (3) 

where pi are the location of the vortex cores. Numeri- 
cally, the minimization of (3) yields an almost triangular 
lattice for a large number of vortex points. 

The aim of this paper is to describe the equivalent of 
(2)- (3) in the case of two-component condensates. Wc 
define gi to be the intra-component coupling strength for 
component i and gi2 to be the inter-component coupling 
strength. For simplicity, wc assume equal masses for the 
atoms in each component and equal trapping potentials, 
but a general case could be handled. The ground state 
of a two-component condensate is given by the infimum 
of 

+512 J i^inv'2p 

under /IV'iP = Ni, /|-02p = ^^2- Wc set gi ^ aig, 
92 = Q^25, 912 = «o5 where g is large, so that e = 
1/y/g is small. We change wave functions to ipi{x,y) = 
y/eui{xy/e,yy/e), i>2{x,y) = y/eu2{xy/e,yy/e). Calling 
a = (ofQ, Q!i, 0:2), the energy wc want to minimize is 



Ea,n{ui,U2) = 




— X r(TOi, Vtti) 



+ '-^\^U2\' + lr'\u2f + ^\u2\' 
— eil X r(TO2, VU2) 
+ ao|uiP|u2p 

(5) 

where (iu, Vw) = iuVu — iuVu. For = 0, the ground 
state is real valued and we denote it by (?7i, 772). It is a 
solution of 

-e^Ar^i + r^ryi -I- 2airil + 2Q:o??2'n = l^-i'Hi (6) 

-e^A772 + r'^m + '2a2r]l + 2aof7?f72 = M2f72- (7) 

The shape of the ground state varies according to a and 
when ttQ — aia2 < 0, can be either two disks or a disk 
and an annulus, as we will see below. 



REDUCTION TO THE CORE EQUATIONS 

We recall that (771, 772) is the ground state for = 
and we consider (ui, 7/2) a ground state of Ea^Q. We call 
(/i:/2) such that ui — 771/1 and 7i2 = 772/2. We expect 



rji to include the slow varying profile and fi to include 
the vortex or peak contribution, so that fi is 1 almost 
everywhere except close to the vortex and peak cores. 
We want to write the energy of (7x1,7/2) as the energy 
of (771,772) plus a rest, which is the energy that we are 
going to study. This follows a trick introduced in [13], 
and used for single Bose Einstein condensates in [9, 10]. 
We multiply (6) by 77i(|/i|2 - 1) and (7) by 772(1/2^ - 1), 
and integrate and add the two equations, which yields 
the identity 

/ y|V7/i|2(|/i|2 - 1) +£\i/iV77i • V/i 
+ ^rWiilM' - 1) + c^ivtm' - 1) + c^ovWiilhl' - 1) 

+ Y|V772 p(|/2p - 1) +e'7;2/2V772 • V/2 
+ ^^'%'(l/2p-l) + «277|(l/2p-l) + aor7|'7?(l/2p-l) = 0. 

(8) 

Note that the Lagrange multiplier term has disappeared 
because Ui and 77^ are normalized similarly. We replace 
(7x1,7x2) by (/1771, /2772) into the energy (5), use the iden- 
tity (8) and find 

Ea,n{ui,U2) = Ea^o{i]i,r]2) + Fa^nifi, f2) where 
i^a,o(/i,/2) = J ^v!\^h\^~eijjnxr{ih,\/h) 

+ la^rJti\h\' 1)' + aovhli^ - l/iHll - IM') 

+ y %'|V/2p - £77^^^ X r(7/2, V/2) + ia2772^|/2p - 1)'. 

(9) 

This splitting of energy does not assume anything about 
the scales of energy: it is an exact identity. We point 
out that as soon as Q;ia2 ~ <Xq > 0, then the quadratic 
form in the energy Fa^n is positive and minimizing Ea^n 
in (7ii,7X2) amounts to minimizing FaS2 in (/i,/2). 



Now we assume that we scale everything close to a 
point p where rj^ ~ pi, rj2 = P2, and /i , /2 can be written 
as functions of p + |?' — p|/e. Then, in the new variable 
f = |r — p\/ e, the functions /i, /2 are a ground state of 

^a.o(/i,/2) = j \pi\Vh?-epinxr[ih,Vh) 

+ \o^lPl{\h? - I? + «0PlP2(l - |/lP)(l - \f2?) 

+ \p2\yf2? - ep2^ X r(7/2, V/2) + \a2pl{\f2? ~ I? 

(10) 
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and solve the system 

- PiA/i - len X rpiV/i + 2a^pl{\h\'' - 
+ 2aoPiP2/i(|/2|'-l) = Ai/i 
- P2A/2 - ien X rp2V/2 + 2a2P2(l/2|' - l)/2 

+ 2aoPlP2/2(|/l|'-l) = A2/2. (11) 

This is exactly the system studied in [4, 14] for a homo- 
geneous condensate. The splitting of energy has allowed 
us to reach a homogeneous system. Assuming a vortex 
in component- 1 and a spike in component- 2, we have 
/i = vi{r)e'^^ and /2 = V2[r). We expect that vi,V2 tend 
to 1 at infinity so that Ai = efi and A2 = 0. This yields 
the following system 



To begin, assume that both components are circular with 
radii i?i and R2 and with i?i < R2. When aia2 — Q^o — 0' 
that is T12 > 0, and e is small, both components are in 
the Thomas-Fermi (TF) regime, and the density profiles 
for r < Ri are 



I |2 1 f "0 2t. 

Zaii 12 \ a2 

I 12 ^ ( 2t 

2a2l 12 V Q^i 



and for Ri < r < R2 are 



A*2 - 
2a2 



(19) 
(20) 

(21) 



r 



^ + 2aipi(wf - l)wi +2aoP2Wi(w2 - 1) =0(12) 



(^^4)' 



with = 0. The chemical potentials /ii and /i2, and 
the radii, Ri and R2, are to be found. In addition we 
have the normalisation condition 



2a2P2(w2 - 1)^'2 + 2aopiV2{vl - 1) = 0.(13) 



From this system, asymptotic expansions can be obtained 
for vi and V2 at infinity: Vi{r) — 1 ^ ^7i/^^ and V2{r) — 
1 ^ 72/f^ for some constants 71 and 72. Equations (12)- 
(13) at infinity imply that Q!2/3272 — olqPi^i and 1 — 
4ai/9i7i + 4q:o/J272 = 0, thus 



\riu\''d = Nk, (22) 

where, for generality, Ni ^ N2. We denote ctk = N^ak 
and q;o = \/NiN20to and get 

/45iri2^ ^'^ 



1 



71 



4/3iairi2 



and 72 



4:p2aia2Ti2 



where 



In particular. 



ri2 



ao7i7iPiP2 



ceia2 



1 -ri2 



(14) 



(15) 



(16) 



Ri 



R2 = 



Ml 



4(52 + (ii(l-ri)) 



1/4 



(23) 
(24) 



/ 4airi2r2 

+(i-r2 



1/2 



[a2 + ai{l - Fi)] 



1/2 



M2 = - [52 + ai(l - Ti)] 



In order to fully analyze the system (12)-(13),we need to 
have information on pi , p2, that is the ground states r/i , 
r]2' in particular, we need to know whether the supports 
of ^1- '?2 arc disks or annuli, and where they reach their 



1/2 



We find from (19) and (14) that 

T2Ri I F27V1 



(25) 
(26) 



Pi = (0) 



2airi 



TraiFi 



(27) 



maximum. 



THOMAS-FERMI PROFILE OF THE GROUND 
STATE 

Wc recall some properties of the solutions of (6)-(7) ob- 
tained in [3]. The following non-dimensional parameters 
are introduced 



and 



while 



71 



16iViaiF2Fi2 



Ti 
F2 



Oo 

a2 



(17) 
(18) 



P2 = tjKO) = — (i-{N2a2 + iViao))'/' - ao{^^^)^/^ 
a2 \ TT Trail 12 

(28) 

and 72 follows from (14). 

Equations (23)- (24) are vahd provided F12/F2 > (to 
ensure that i?i and fii are real) and aiFi < ci2F2 (to 
ensure that R2 > Ri). If instead the initial assumption 
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on the size of the radii was taken to be i?i > R2, then 
the appropriate expressions would also be given by Eq.'s 
(23)- (24), however with the indices 1 and 2 alternated. 
In this case, the conditions would be fiiFi > 0:2^2 and 

ri2/ri > 0. 

Returning now to R2 > one must ensure that 
|?72p > for all r < i?i to have a disk rather than an an- 
nulus. Suppose that there is a point at the origin, where 
|r;2p = 0. Then, from Eq. (20), 

aifJ.2 
Ml 

Niai 1 I alNf 4N2aia2 , . 

2{Ni + N2) 2 Y (iVi + N2y ivTHv 

The existence of some at which the density in 
component-2 hits zero at the origin is the indication 
of a spatial separation of the components. Notice that 
this critical value for ri2 is independent of il. In the 
spatial separation regime, component- 1 is circular while 
componcnt-2 is annular, provided R2 > Ri. It is not 
possible for an annulus to develop in component- 1 if 
i?2 > Ri', this can be seen by writing down the TF den- 
sity expressions for an annular component- 1 and a circu- 
lar component-2 in which the chemical potentials become 
multi-valued. Similarly, an annulus can only develop in 
component- 1 if Ri > i?2- Thus under the assumption 
R2 > Ri, an annulus can only develop in component-2 
and the condition to have two disks is thus ao < cio- 



componcnt-2, then they minimize the point energy 
- Tre^ pi log -Pj\- TTE^ p2 log \q^ - q, \ 

+ "h'"l°S^"2^+^"^0^"^^"' 

^ ^ i 

+ -(-e'noge\^+snp2)Yh\^ 

Estimate of the kinetic energy term 

Let us call pi the vortices in component- 1, and 
Qj in componcnt-2. Then the kinetic energy term 
(1/2) J 77ie^|V/ip provides a leading order term due to 
the kinetic energy of the phase (which behaves locally 
like 1/r outside a disk of radius e around each vortex) , 
which is 

Tre^ Y PTF, 1 (Pi ) I log £ I - 7re^ ^ ptf, 1 {Pz ) log \pi - Pj \ 

i i=jtj 

(34) 

with a similar term for component-2, where Pi is replaced 
by qy 

Rotation term 



We call Xi{r) the primitive of — ?'Ptf.i(^) which van- 
ishes at i?i and X2{r) the primitive of —rpTF,2{f) which 
vanishes at R2. Then for r < Ri 



VORTEX INTERACTION 



-iRl~ry 



(35) 



Let us call pt_f,i, Ptf,2 the Thomas-Fcrmi limits of 
|77i|2 and |r,2|2 given by '(19)-(20)-(21) and (23), (24), 
(25), (??), (26) in the case of two disks. Then 



PTFS 



2aiFi2 



[Rl 



(30) 



PTF,2 = 7r^^{Rl~r^) + 7^{Rl^Rl) ifr<R, (31) 
2a2i 12 2a2 



and {Rl-r^) if r > 

2a2 



(32) 



We want to estimate the various terms in the energy Fa^n 
as in [9, 10] and we are going to show that, if pi arc 
the vortices for component- 1 and qj are the vortices for 



X2{r) = ^^^{Rl-r'r + ^{Rl-Ri){Rl+Rl-2r'). 
oa2^ 12 oa2 

(36) 

Thus the rotation term —eft j rjl y. r(i/i,V/i) is well 
approximated by —erj / VXi x (i/i, V/i). An integration 
by parts around each vortex yields 

-2TTenY,Xi{pi) (37) 

i 

with a similar contribution for component-2. 

First vortices 

The leading order approximation of the kinetic and 
rotation energy yields (assuming vortices at points pi for 



5 



component-1 and qj for componcnt-2): 

7re^ X] PTF.i{Pi)\ log e| + tte^ ^ ptf,2 (gj ) | log £| 

i 3 

-2T:eVt'Y^Xi{pi)^2T:en'Y^X2{qj). (38) 



The energy is minimized by a configuration such that 
PTF,ilXi or ptf^iI^i reaches its minimum. We find 
that 



PTF,i{r) 4 



(39) 



1.^2 _ 2^ , 1 (i?i-flf)(fli-r^) 



(40) 

This imphes that above a critical value fi'^, vortices be- 
come energetically favorable in the system and VL'^ is given 
from (38) by 



-e|log£|mm— — - 

Z i.r A,- 



(41) 



For a harmonic potential, and in the case of two disks, 



the minimum of occurs at the origin. Furthermore, 



since 



where 



Ptf.2 _ Ptf.i 



-D{r) 



(42) 



A{Rl-Rl){Rl-r-) 
^^'>~ {Rj r2) ^ 



if^(i?f - + {Rl - Rim + Rl- 2r2) 

(43) 

and given the signs of the parameters, we see that D(0) > 
so that we always have ptf,i/Xi > ptf,2/X2- 



The first vortex is thus preferred in componcnt-2 (i.e. 
the component with larger support) and occurs at the 
origin with the critical velocity given by 
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FIG. 1: The critical velocity for creation of the first vortex 
plotted analytically from Eq. (44) [solid line] and numerically 
[dotted line] as a function of for two parameter sets: (a) 
e = 0.0352, ai = 0.97, 02 = 1.03 and (b) e = 0.0358, ai = 

Q!2 = 1. 



Since miui^r PTF,i/ Xi = 4:/Rl, we have 



= ejlogel 



ttFi 



aiF 



12 



(45) 



We have plotted 51^ — F12 curves for two cases in 
Fig. 1, the first with distinct intracomponent coupling 
strengths and the second with equal intracomponent cou- 
pling strengths (where Eq. (44) reduces to Eq. (45)) and 
compared then to the numerical results of [3] (these pa- 
rameter sets correspond to sets 'ESI' and 'ES3' respec- 
tively from [3]). We find good agreement between the 
two results. 



17" 



a2Fi 



logej fFi2V"2 + ai(l -Ti) 



(l-FOv/aiFsFi 



(44) 



Note that this expression gives H'^ = when ao = ao 
(provided ai 7^ a2 otherwise W reduces to a non-zero 
constant) . For some computations in the rest of the pa- 
per, we can assume TViai = N20i2 so that i?i = R2 and 
we have a lattice of peaks and vortices close to the origin. 



Energy expansion 



We now have to go further into the energy expansion 
to estimate the interaction energy. We assume that the 
vortices appear close to the origin. Then (34) and (37) 
can be expanded around the origin, using (30)-(32) and 
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(35)- (36), which yields 



^ Pl log \p^ ~Pj\~ TTE^ ^ p2 log |g, - qj I 

r2 



Tre —el log el 



2airi2 



(46) 
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Interaction energy 



We find from (9) that the interaction energy is 

aopiP2 [ (i-iwin(i-iw2n. 



Near a vortex-peak, this reduces to 

1 1 



aoPiP27i72 



^(1,0) ^^0,1) 



where we take the notations of [4]: rji q) is the local 
distance to the vortex in componcnt-l, and 7^(0.1) is the 
distance to the next peak in component- 1, or cquivalently 
to the vortex in component-2. From (16), we find that 
the coefficient in front of the integral is equal to .It^-M ■ 

^ ^ 16ail ^2 

The computations in [4] allow us to estimate the integral 
term and we find for the interaction term 

i-ri2 1 41, I 1 

16rf2 ai |p, - gjl^ 

This is for a vortex in vi. Of course, if the vortex is in 
V2, it would be different by a factor l/a2- 

The interaction energy is thus 



i~ri2 
i6rf2 



1 

ai 



a2 



(47) 



Together with (46), this leads to (33). 



NUMERICAL SIMULATION OF THE 
RENORMALISED ENERGY 



FIG. 2: The ground state of the renormalised energy (50) 
calculated for (a) c^i — 0.05 and (b) c^i = 0.3 with the number 
of points in each component equal to = 46. Component- 1 
is shown in the left columns and component-2 in the right 
columns. 
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FIG. 3: The critical value of Q^" as a function of calculated 
analytically by Eq. (55) (solid line) and numerically (dashed 
line) for the parameters e = 0.0358, ai = 02 = 1. 



We now want to find the ground state of (33) when the 
radii of both components are equal. Under this condition 
we have pi — P2, ^1 = ^2 and ai ~ a2- This allows us to 
perform a rescaling that leaves the renormalised energy 
only dependent on a single parameter. We write f2 = wfi"^ 



for i}" defined in (45). Then 

(48) 
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which implies that we can rescale the pi's and gj's as 
Pi = ^pi and qj = jq^ with 




so that the new energy is 



We plot this form of fi** as a function of T12 in Fig. 
3 where we have taken s = 0.0358 and ai — a2 = 1 
(note that this parameter set corresponds to set 'ES3' in 
[3]). This provides good agreement with simulations of 
the full Gross Pitaevskii energy and confirms that our 
point energy (33) well describes the system. 



with 



7r(l-ri2) 

4air?2 



, (2^-l)e2|loge|^ 



(50) 



(51) 



We simulate this renormalised energy (50) using a con- 
jugate gradient method varying the parameter c^^ and the 
number of vortex points. For a fixed number of lattice 
points A^. when increasing c^j, the ground state lattice 
goes from triangular to square at a critical c^^ . Note that 
when N gets large, c*^ no longer depends on N . In Fig. 
2 we plot the ground state for two values of Cu which give 
a triangular [c^ ~ 0.05) and a square lattice [cu = 0.3) 
when TV = 46. 

For each iV, we can calculate the critical c^^ and com- 
pare this with the simulations on the full GP equations, 
as performed in [3] (they provide the appropriate value 
of N). From (51), we thus find the critical value of f2 for 
which the lattice goes from triangular to square: 



1 



ttFq 



2 ^(l-ri2)e2|log£n V ^^1^12 



logel (52) 



It turns out that when fi gets close to 1, the condensate 
expands and one has to include in the TF profile a term 
(1 — ri2)r2 instead of just r^. The radii i?i and R2 vary 
like (1 — $12)1/4^ This changes from (51) to 



7r(i-ri2) / 2w- vT^TP 



I logel' (53) 



so that if we define 



/3t 



2 ^(i-ri2)|iog£r / V 5iri 



then = Pts \J\ - (Sl^' f, which yields 



I logel (54) 



(55) 



DISK PLUS ANNULUS 

Equation (29) introduces the critical aQ (or equiva- 
Icntly ri2) for which the two disk state is no longer a 
solution. In this section we consider ag > ao so that 
the system is given by a disk in component- 1 and an an- 
nulus in component-2. This requires ai to be different 
to a2, because otherwise the conditions ao > ao and 
ap — aia2 < are not consistent. Defining the inner 
and outer radii of componcnt-2 to be at r = i?2_ and 
r = i?24. such that i?2_ < i?i < i?2+ , the appropriate 
density profiles are 



PTF.iir) = 1 771 1 2 
for < 7' < i?2_ , 



Ml 



2ai 



(56) 



PTF.iir) = I771I 



' - ^ ^ii-%2-r2r2)(57) 



2airi2 V a2 



PTF,2{r) = \ri2\ 



1 



ao 



2a2Ti2 \ ai 



M2 - —Ml - r'Ti (58) 



for i?2_ < r < Ri, and 



t ^ I ,2 P2-r^ 

2a2 



(59) 



for i?i < ?■ < i?2+ . 

The normalisation condition (22) gives the following 
expressions for the radii 



R2 




(60) 



Fi 



Ri = \iR-\i-Y-^a~^2rs (61) 



and 



i?2+ = yR + v-^s 

Pi = R 

P2 = R+ ^-ViT2S 



(62) 

(63) 
(64) 



for the chemical potentials. Here we have introduced the 
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parameters 



R 
S 



= 2\ 



= 2, 



/ai(l + |f(l-r2)2) 



(65) 
(66) 



A number of conditions can immediately be found from 
these expressions. In order for i?2_ < it must be that 
> [The inequahties i?2+ > Ri and i?2^ > i?2_ are 
automatically valid]. Secondly, all the expressions (ex- 
cept that for ^i) require rir2 < 0. This is equivalent to 
Q!2 > cti (under the assumption that the annulus devel- 
ops in component-2). Thus the range of ri2 for which 
component- 1 is circular and component-2 is annular is 
given by < ri2 < ri2, where ri2 = ri2 evaluated at 
ao = Q^o given by (29). 

When ao ~ we have shown in Eq. (44) that the 
critical velocity for creation of the first vortex (that oc- 
curs at the origin in component-2) is identically zero. 
This is precisely the point at which the annulus develops 
{pTF,2{'r = 0) = with i?2_ = 0). In the region ag > aoi 
the first appearance of a topological defect is the devel- 
opment of a giant vortex in component-2 and this can be 
analyzed fully. We need to introduce 



As 



02 
0-2 



PTF,2{s) 



R2 



4a2ri 



rii?2_ log 



ds 



Rl 



ri2i?^^ log 



Rl 







Then the critical velocity for nucleation of the giant vor- 
tex is determined by Vlg^ = £^2 and the circulation of 
this giant vortex is the integer part of f2/(eA2). Next, we 
want to determine vortices in the bulk and apply (38) to 
the case of the disk plus annulus. Thus we need to define 
Xi(r) to take into account the giant vortex. 

We define ^\{r) the primitive of ~rpTF,iif) which 
vanishes at Ri giving 



Xiir) 



1 



2\2 



[Ri-r^) 



8aia2ri2 



Rim 



R'. 



2r2) 



(68) 



for < r < i?2_ , and is the same as Eq. (35) in the 
region i?2_ < r < Ri. In order to account for the de- 
velopment of vortices beyond the giant vortex, we must 
recall that additional to the term —efi / '^Ptf,2(«/2, V/2), 
there will be a kinetic energy term coupling the giant vor- 
tex and the vortex cores, namely e^d / Ptf,2/'^(«/2, V/2), 
where d is the degree of the giant vortex, of order 
il/{eA2). Therefore, we define X2{r) to be the primi- 



tive of —rpTF,2{r) + (1/A2)/Otf.2 ('')/?' which vanishes at 
both i?2_ = 0' and R2. = 0, 



X2{r) 



R2. 



spTF,2{s)ds- — 
A2 



1 f'^'+ PtfAs) 



ds (69) 



where A2 is given by (67). Completing the integrals gives 
1 ^^^ aor2 ^^2 _ ^2^2 



1 



4a2ri2A2 



8Q;ia2ri2 
TiRl log 



Rl 



(70) 



ri2i?^^ log ' 



Hi 



- TiRj 



for i?2_ < 7 < Rl and 
X2{r)=^{Rl~r'r 



8a2 + 
1 



4a2A2 



r'-Rl 



(71) 



for Rl < r < i?2+. As before we write down the ratios 
PTF,i/Xi: 

PTF,i 4 _ / 4riao {RL ' RDiRL - r') \ 

Xi ^'■^ (i??-r2) I^ri2a2 (Rj-r^) J 

X '{Rl^r'f 



ri2Q!2 

for < r < i?2^. 



{RL-Rl){Rl + Rl_-2r') 



(72) 



Ptf.i 

PTF,2 
X2 



(r) 
(r) 



(i?f-r2)' 

i 12 



(73) 



(RL-r^) 



2,2 Ct0^2 ir,2 2 



aiEi 



{Ri-r^ 



^ 'rii?Liog^'^^ 



A2ri2 
32 



Rl 



+T12RL log . ^2 



ri(r2-i?L) (74) 



for i?2_ < r < Rl and 

PTF.,2 , 



X2 



4{Rl-r') 



(Rl - r^r h Rl. i«g ^ 



r2 - i?2 
(75) 
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FIG. 4: (a) Functions ptf,i/X\ (solid lines) and PTF.2/X2 
(dashed lines) given by Eq's (72)-(75) for = 0.0202 (black 
lines) and ri2 = 0.015 (blue lines), (b) The critical velocity 
for creation of the first vortex plotted from Eq. (41) [black 
line] and numerically [blue line] as a function of ri2. The solid 
part corresponds to a vortex in component 1 at the center and 
the dashed part in component 2 in the interface region. The 
parameters are e = 0.0352, cti = 0.97 and 02 = 1.03. Note 
that ri2 = 0.02 



poncnt 2. It is not possible to find an analytic expression 
for this minimum and so we must resort to plotting the 
function and comparing this to m.in{pTF.i/ Xi). In Fig. 
4, we plot the functions ptf,i/Xi and ptf,iI^i in two 
different cases where min(pTF,i/-'^i) < "^^"^{ptf,-!! X-i) 
and mm{pTF,\l Xx) > mm{pTF,2/ X2)- We also plot the 
critical Q defined in Eq. (41) as a function of ri2. We 
see that until some critical value of ri2, the vortex first 
nucleates at the origin in component- 1, while above in 
some small region in before reaching ri2, the vor- 
tex first nucleates in the interface region in the second 
component. Then, increasing fJ, leads to a vortex lattice 
close to this point of nucleation. Recall that the giant 
vortex is always present. 



Conclusion 

From the energy of a rotating two component conden- 
sate, wc have derived a reduced energy (33) governing the 
location of peaks and vortices. We have found that the 
ground state of this reduced energy yields a square lat- 
tice of vortices in regimes consistent with the ones found 
numerically for the full Gross Pitacvskii energy. Wc have 
analyzed in detail the geometry of the ground state (two 
disks or disk and annulus) and derived formula of the 
critical velocity for the appearance of the first vortex. 
We can determine in which component the first vortex 
appears and the shape of the lattice. 
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. / Ptf.i 



a2r2fl.i+aorifl.2_ 



for 
for 



< r < i?2 
i?2_ <r < 



Ri 



(76) 



Using the expressions above for the radii, we see that 



a2T2Rl + aoTiRl_ 2 In , . I " , . 

4 p4 I r p4 = — W— (Qi«2 +aia2) (77) 

Q!2l -I- aol i/t2_ ai V "2 

which we note is independent of ag. If we compare the 
values of Eq. (76), we see that they are equal when 
ao = ao (which we can consider a degenerate case) and 
otherwise min(yOTF,i/^i) over the whole space is simply 
given by the expression in Eq. (77). This implies that the 
vortex in component- 1 is preferred at the origin rather 
than in the interface region of the two components. 

To complete the analysis we must also find the 
vam(^PTF,2l X-2) over the whole space to determine 
whether vortices appear first in component 1 or com- 
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